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The degree to which a pure quantum state is entangled can be characterized by the distance or 
angle to the nearest unentangled state. This geometric measure of entanglement, already present in 
a number of settings [A. Shimony, Ann. NY. Acad. Sci. 755, p. 675 (1995) and H. Barnum and N. 
Linden, J. Phys. A: Math. Gen. 34, p. 6787 (2001)], is explored for bipartite and multipartite pure 
and mixed states. It is determined for arbitrary two-qubit mixed states and for generalized Werner 
and isotropic states, and is also applied to certain multipartite mixed states. 
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Introduction: Only recently, after more than half a cen- 
tury of existence, has the notion of entanglement become 
recognized as central to quantum information process- 
ings As a result, the task of characterizing and quan- 
tifying entanglement has emerged as one of the prominent 
themes of quantum information theory. There have been 
many achievements in this direction, primarily in the set- 
ting of bipartite systems 0. Among these, one high- 
light is Wootters' formula H for the entanglement of for- 
mation for two-qubit mixed states; others include corre- 
sponding results for highly symmetrical states of higher- 
dimensional systems 0]. The issue of entanglement 
for multipartite states poses an even greater challenge, 
and there have been correspondingly fewer achievements: 
notable examples include applications of the relative en- 
tropy j^, negativity and Schmidt measure [p^ . 

In this Letter, we present an attempt to face this chal- 
lenge by developing and investigating a certain geometric 
measure of entanglement, first introduced by Shimony |^ 
in the setting of bipartite pure states and generalized 
to the multipartite setting (via projection operators of 
various ranks) by Barnum and Linden We begin 
by examining this geometric measure in pure-state set- 
tings, and then extend it to mixed states, showing that 
it satisfies certain criteria for good entanglement mea- 
sures. We demonstrate that this measure is no harder to 
compute than the entanglement of formation E-p , and ex- 
emplify this fact by giving formulas corresponding to 
for (i) arbitrary two-qubit mixed, (ii) generalized Werner, 
and (iii) isotropic states. We conclude by applying the 
geometric entanglement measure to certain families of 
multipartite mixed states, for which we provide a prac- 
tical method for computing entanglement, and illustrate 
this method via two examples. 

It is not our aim to cast aspersions on exisiting ap- 
proaches to entanglement; rather we simply wish to add 
one further aspect to the discussion. 

Basic geometric ideas; application to pure states: We 
begin with an examination of entangled pure states, 
and how one might quantify their entanglement by 
making use of simple ideas of Hilbert space geome- 



try. Let us start by developing a quite general for- 
mulation, appropriate for multipartite systems compris- 
ing n parts, each of which can have a distinct Hilbert 
space. Consider a general n-partite pure state |?/') — 

I (1) (2) (n), r\ 

Z^pi - p„ XpiP2 - P„|epi • • •e^„'). One can envisage a 
geometric definition of its entanglement content via the 
distance d = miuj^^ || — \<t')\\ between and the near- 
est separable state |0) (or equivalently the angle between 
them). Here, |0) = ^f^il^*-*-*) is an arbitrary separable 
(i.e., Hartree) n-partite pure state, the index i = 1 . . .n 

labels the parts, and = J2pi '^Pi \^Pi)- Naturally, 

the more entangled a state is, the further away it will 
be from its best unentangled approximant and the wider 
will be the angle between them. 

To actually find the nearest separable state, it is conve- 
nient to minimize, instead of d, the quantity || — 
subject to the constraint {4>\(j)) = 1. In fact, in solv- 
ing the resulting stationarity condition one may restrict 
one's attention to the subset of solutions \(j>) that obey 
the further condition that for each factor If^''-*) one has 
= 1. Thus, one arrives at the nonlinear eigen- 
problem for the stationary \(f)): 

5Z ^PlP2---Pn'^Pl ' ' ' ^Pi ' ' ' ^pj ~ ^^Pi ' (-'-^) 

Pl---Pi---Pn 

y ] XpiP2---pn'^pi ' ' ' Cpi • • • CpJ — Acp.-* , (lb) 

Pl---Pi---Pn 

where the eigenvalue A is associated with the Lagrange 
multiplier enforcing the constraint — 1, and ^ de- 

notes exclusion. In basis-independent form, Eqs. (^ read 

I |0(^)))=A(0«|, ( ® (0(^)|)|^)=A|0«).(2) 

From Eqs. (|l|) or (|) one readily sees that the eigenvalues 
A are real, in [—1, 1], and independent of the choice of 
the local basis {|ep'^)}. Hence, the spectrum A can be 
interpreted as the cosine of the angle between and 
|(/)); the largest, Amax, which we call the entanglement 
eigenvalue, corresponds to the closest separable state. 
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Although, in determining the closest separable state, 
we have used the squared distance between the states, 
there are alternative (basis-independent) candidates for 
entanglement measures: the distance, the sine, or the 
sine squared of the angle 9 between them (with cos 6* = 
'Re {tp\(j))) . We shall adopt iSgin^ = 1 — ^max 
tanglement measure because, as we shall see, when gen- 
eralizing i?sin^ to mixed states we have been able to show 
that it satisfies a set of criteria for entanglement mea- 
sures. We remark that determining the entanglement 
of IV') is equivalent to finding the Hartree approxima- 
tion to the ground-state of the auxiliary Hamiltonian 

In bipartite applications, the eigenproblem (|l|) is in fact 
linear, and solving it is actually equivalent to finding the 
Schmidt decomposition [Q. Moreover, the entanglement 
eigenvalue is equal to the square of the maximal Schmidt 
coefficient. By constrast, for the case of three or more 
parts, the eigenproblem is a nonlinear one, for which one 
can, in general, only address the problem directly, i.e., by 
determining the eigenvalues and eigenvectors simultane- 
ously, presumably numerically. Yet, as we shall illustrate 
shortly, there do exist certain types of states whose en- 
tanglement eigenvalues can be determined analytically. 

Illustrative examples: Suppose we are already in posses- 
sion of the Schmidt decompostion of some two-qubit pure 
state: — ^/p |00) + ^Jl — p\ll). Then we can read off 
the entanglement eigenvalue: A^ax — uiaxjy^, y'l — p}. 
Now, recall that the concurrence C for this state is 
2\/p(l — p)- Hence, one has 

ALx = Hi + V^~C^ : (3) 

which holds for arbitrary two-qubit pure states. 

The possession of symmetry by a state can allevi- 
ate the difficulty associated with solving the nonlin- 
ear eigenproblem. To see this, consider a state \tp) — 

Spi...p„ XpiP2-P„l4i^42^ • • • 4"^> that obeys the symme- 
try that the nonzero amplitudes x are invariant under 
permutations. What we mean by this is that, regardless 
of the dimensions of the factor Hilbert spaces, the am- 
plitudes are only nonzero when the indices take on the 
first V values (or can be arranged to do so by appropri- 
ate relabeling of the basis in each factor) and, moreover, 
that these amplitudes are invariant under permutations 
of the parties, i.e., Xo-io-2 - <t„ = XpiP2 - pnJ where the ct's 
are any permutation of the p's. (This symmetry may 
be obscured by arbitrary local unitary transformations.) 
For such states, it seems reasonable to anticipate that 
the closest Hartree approximant retains this permutation 
symmetry. Assuming this to be the case — and numerical 
experiments of ours support this assumption — in the task 
of determining the entanglement eigenvalue one can start 
with the Ansatz that the closest separable state has the 
form 10) = Cj\e^^'^)) , i.e., is expressed in terms 

of copies of a single factor state, for which c^*^ = Cj. To 
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FIG. 1: Entanglement of the pure state -/s | W) + \/l — s | W) 
versus s. This also turns out to be the entanglement curve 
for the mixed state s |W){W| + (1 - s)|W){W|. 

obtain the entanglement eigenvalue it is thus only nec- 
essary to maximize Re with respect to {cjj^^j^, a 
simpler task than maximization over the X^ILi ampli- 
tudes of a generic product state. 

To illustrate this, we consider several examples involv- 
ing permutation-invariant states, restricting attention to 
V = 2. The most natural realizations are n-qubit sys- 
tems. One can classify these symmetric states, as follows: 

\S{nM^\[^^^ E I0_01_l)- (4) 

permutations ^. n—k 

As the amplitudes are all positive, one can assume that 
the closest Hartree state is of the form \(f)) = (v^|0) + 

^1 — p jl))*"", for which the maximal overlap (w.r.t. p) 
gives the entanglement eigenvalue for |S(n, fc)): 

For fixed n, the minimum Amax (and hence the maximum 
entanglement) among the |S(n, fc))'s occurs for k = n/2 
(for n even) and k = (n ± l)/2 (for n odd). In fact, for 
fixed n the general permutation-invariant state can be 
expressed as ak\S{n, k)) with J2k lo^fcP — 1- The en- 
tanglement of such states can be addressed via the strat- 
egy that we have been discussing, i.e., via the maximiza- 
tion of a function of (at most) three real parameters. 
The simplest example is provided by the nGHZ state: 
InGHZ) = (|S(n,0)) + |S(n,7i)))/V2. It is easy to show 
that (for all n) Amax("GHZ) ^ 1/V2 and £'si„2 = 1/2. 

We now focus our attention on three-qubit settings. Of 
these, the states |5(3,0)) = |000) and |5(3,3)) = |111) 
are not entangled and are, respectively, the components 
ofthe the 3-GHZ state: |GHZ) = (|000)-F|lll))/v^. The 
states |S(3,2)) and |S(3, 1)), denoted |W) = |S(3,2)) = 
(|001) + |010) + |100))/y3and |W) = |S(3, 1)) = (|110> + 
|101) -I- 1011))/\/3, are equally entangled, having Amax = 
2/3 and £'si„2 = 5/9. ^ 

Next, consider a superposition of the W and W states: 
|WW(s,0)) = Vi|W) + x/r^e^-^lW). It is easy to see 
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FIG. 2: Entanglement of \WG{s,(f))) versus s. The upper 
curve is for (j) — n whereas the lower one is for = 0. Dots 
represent states with randomly generated s and 0. 

that its entanglement is independent of 0: the transfor- 
mation { |0), |1) } ^ {|0),e-*'^|l)} induces |WW(s,0)) ^ 

e-"^|WW(s, 0)). To calculate A^ax, assume that the 
separable state is (cos6'|0) -I- sin6'|l))'^'^, and maximize 
its overlap with |WW(s,0)). Then we find that the 
tangent t = tan 6* is the particular root of the poly- 
nomial Vl - si^ -f - 2Vl - st - = that 
lies in the range t g [-^1/2, \/2]. Via 9{s), Amax (and 
Esin^ 1 — ^max) <^an bc cxprcsscd as 

A,„ax(s) = i(Vscos6'(s)-|-VT^sin6l(s)) sin26l(s). (6) 

In Fig. 1^, we show E^in^ (|WW(s, (f>))) vs. s. In fact, A,„ax 
of the more general superposition 



\SSn-Mk, {r, cf>)) = V7 |S(n, h)) + yi^e*"* \S{n, fca)) (7) 

(fci ^ turns out to be independent of 0, as in the 
case of |WW(s,^)), and can be computed in the same 
way. We note that although the curve in Fig. |l| is convex, 
convexity does not hold uniformly over fci and ^2. 

For our last pure-state example, we consider superpo- 
si tions of W and GHZ states: |WG(s,0)) = ^/s\SN) + 
Vl — s e*'^|GHZ). For these, the phase <j) cannot bc 
"gauged" away and, hence, -Esin^ depends on (/). In Fig. |^ 
we show i?sin2 vs. s at = and tt (bounding curves), as 
well as i?sin2 for randomly generated values of s G [0, 1] 
and (j) S [0,27r] (dots). It is interesting to observe that 
the 'tt' state has higher entanglement than the '0' does. 
As the numerical results suggest, the (0-parametrized) 
i?gi„2 vs. s curves of the states |WG(s,(/>)) lie between 
the 'tt' and '0' curves. 

Extension to mixed states: The extension to mixed states 
p can be made via the use of the convex roof (or hull) con- 
struction [indicated by "co"], as was done for the entan- 
glement of formation (see, e.g., Ref. |]). The essence is a 
minimization over all decompositions p = 'Y^^Pi 
into pure states, i.e., 

E{p) = (coEp^^c){p) = min V p., ^;pure(|V'j))- (8) 

Now, any good entanglement measure E should, at least, 
satisfy the following criteria (c.f. Refs. 0, E2I Esl): 



CI. (a) E{p)>{)] (b) E[p) = if p is not entangled. 
C2. Local unitary transformations do not change E. 
C3. Local operations and classical communication 
(LOCC) (as well as post-selection, if one wishes) do not 
increase the expectation value of E. 
C4. Entanglement is convex under the discarding of in- 
formation, i.e., J^iPi^iPi) > E{J2iPiPi)- 
The issue of the desirability of additional features, such as 
continuity and additivity, requires further investigation, 
but C1-C4 are regarded as the minimal set, if one is to 
guarantee that one has an entanglement monotone [ p^ . 

Does the geometric measure of entanglement obey Cl- 
4? From the definition (||) it is evident that CI and C2 
are satisfied provided that -Epuro satisfies them, as it does 
for -Epuro being any function of Amax consistent with CI. 
It is straightforward to check that C4 holds, by the con- 
vex hull construction. The consideration of C3 seems 
to be more delicate. The reason is that our analysis of 
whether or not it holds depends on the explicit form of 
^'puro- For C3 to hold, it is sufficient to show that the 
average entanglement is nonincreasing under any trace- 
preserving, unilocal operation: p — > J^k^^pV^ , where 

the Kraus operator has the form = 1 ® • • • 1 ® V^^*' <Si 
1 • • • (g) 1 and J2k ^k^k = 1- Furthermore, it suffices to 
show that C3 holds for the case of a pure initial state, 
i.e., p = \'ip){ip\. We now prove that for the particu- 
lar (and by no means unnatural) choice -Epurc = ^-sin^, 
C3 holds. To be precise, for any quantum operation on 
a pure initial state, i.e., |V')(V'I J2k^k\4'){'4'\^k ^ we 
aim to show that J2kPk E^in^ {yk\i>) /^/Pk) < £'si„2 (|f/')), 
where pk = Tr Vfc|?/')('0| = {i/j\V^Vk\ilj) , regardless of 
whether the operation {Vk} is state-to-state or state- 
to-ensemble. Let us respectively denote by A and A^ 
the entanglement eigenvalues corresponding to and 
the (normalized) pure state Vfe|'!/')/y^fe . Then our task 
is to show that J^kPk^k — '^^ which the left 

hand side is, by the definition of A^, equivalent to 

EkPk.TojmiVkW/vp^r = EkC^imiVkm'- 

Without loss of generality, we may assume that it is the 
first party who performs the operation. Recall that the 
condition (||) for the closest separable state \4>) = \a)ii^ 
17)2.. .„ can be recast as 2---n(7|''/')i---n=A|a)i. Then, 
by making the specific choice {^k\ — ((^l V^^^^VV^fc) 



(7I, where qk = ((5|Vj 



(l)tT/(l)|;^.\ 



sought result Y^kPk^l = E 

i2Y^ ^/;T,n/(l)tT/(l)|;;,\ / /77^^2 



'V^'''\a), we have the 



1(61^1^)11 



> 



^'Ekmvryn^)/v^kr 



max 

k ikeD 

= A2. We note that a dif- 
ferent approach to establishing this result has been used 
by Barnum and Linden 

Before moving on to the terra incognita of mixed mul- 
tipartite entanglement, we test the geometric approach in 
the setting of mixed bipartite states, by computing ii^sin^ 
for three classes of states for which Ep is known. 
Arbitrary two-qubit mixed states: For these we show that 



E ■ ' 



■ip) 



C{P? 



(9) 
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FIG. 3: Entanglement curve for the mixed state p7-2,5{r) (full 
line) constructed as the convex hull of the curve for the pure 
state \SS7;2,5{r,(j))) (dashed in the middle; full at the edges). 

where C{p) is the Wootters concurrence of the state p. 
Recall that in his proof of the formula for Ep, Wootters 
showed that there exists an optimal decomposition p — 
X^i^j in which every lipi) has the concurrence 

of p itself. (More explicitly, every \tpi) has the identical 
concurrence, that concurrence being the infimum over all 
decompositions.) By using Eq. (||) one can, via Eq. (||), 
relate i^sin^ to C for any two-qubit pure states. As i^sin^ 
is a monotonically increasing function of C G [0, 1], the 
optimal decomposition for -EginS is identical to that for 
the entanglement of formation Ep. Thus, we see that 
Eq. (^) holds for any two-qubit mixed state. The fact that 
i?si„2 is related to E-p via the concurrence C is inevitable 
for two-qubit systems, as both are determined by the one 
independent Schmidt coefficient. 

Generalized Werner states: Any state pw of a C'' (g) 
system is called a generalized Werner state if it is in- 
variant under Pi : p ^ J dU{U U)p{U'^ ® [/^), where 
U is any element of the unitary group hi{d) and dU is 
the corresponding normalized Haar measure. Such states 
can be expressed as a linear combination of two opera- 
tors: the identity 1, and the swap F = 'Yliij i-S-, 
= al + &F, where a and h are real parameters re- 
lated via the constraint Trpw = 1- This one-parameter 
family of states can be conveniently expressed in terms 
of the single parameter / = Tr(pvvF). By employing the 
technique by developed by VoUbrecht and Werner |6| [as 
applied to i?F(/3w)] to i?sin='j one arrives at the geometric 
entanglement function for Werner states: 



and zero otherwise. 

Isotropic states: These are states invariant under P2 : 
p ^ J dU {U ®U*)p {W ®U*^), and can be expressed as 



Piso{F) 

where |<i>^ 



1 - F 



d2 



- (i -[$+>($+ 1) + F|ci>+)(ci>+ 1, (11) 



^ T^Etil") and F e [0,1]. For F e 
[0, 1/d], this state is known to be separable [Q. By us- 
ing the Vollbrecht- Werner technique and following argu- 
ments similar to those of Terhal and Vollbrecht ap- 



plied to i?F(piso), one arrives (for F > 1/d) at 

Esin^- (pUF)) = 1 - ^{VF+Vi^-P)id-l)y. (12) 

Mixtures of multipartite symmetric states: As a final ex- 
ample we consider mixed states of the form (fci 7^ ^2) 

Pn-MkJ(>') = T |S(7i, A;i))(S(7i, A;i)|-('(l-r)|S(n, ^2)) (S(n, A:2)|- 

From the independence of £^si„2 (|SS„;fejfe2 ('^i 0))) on 0, 
one can show that Egi^2 [pn-Mk^i''')) vs. r can be con- 
structed from the convex hull of the entanglement func- 
tion of |SS„;fejfe2 (^i 0)) vs. r. An example is shown in 
Fig. H. If the dependence of E^^^-^2 on r is already con- 
vex for the pure state, its mixed-state counterpart has 
precisely the same dependence. Figure |l|, for which 
(n, fci, /c2) — (3, 1, 2), exemplifies such behavior. 

Concluding remarks: We have considered a general, geo- 
metrically motivated measure of entanglement, applica- 
ble to pure and mixed quantum states involving arbitrary 
numbers and structures of parties. We have illustrated 
this measure via several examples. In bipartite settings, 
this approach provides an — in general, inequivalent — 
alternative to the entanglement of formation; it is, more- 
over, naturally extendable to multipartite settings. 
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